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Abstract 

In models of modified gravity, extra degrees of freedom usually appear. 
They must be removed from the spectrum because they may indicate the 
presence of instabilities and because otherwise the model might not agree 
with observation. In the present letter, we will discuss a model that modifies 
General Relativity through the addition of a Lorentz-violating potential-like 
term for the metric. No additional propagating modes and no classical insta- 
bilities are present. The model departs from GR only in the extreme infrared 
region, and the cosmological evolution contains a bounce when the size of the 
Universe is of the same order of the added deformation. 



1 Introduction 



The gravitational interaction is well described by General Relativity within a vast 
range of distance scales: from solar system to galaxy cluster sizes. Still questions 
arise on its applicability at microscopic, where quantum effects are no longer negli- 
gible, and super-horizon distances. At the far end of the length scale, in the deep 
infrared region, it is not obvious whether gravity should behave the same as it does 
at much shorter distances. Rather, the discovery of (recent) cosmic acceleration 
from supernovae data opens up the possibility that gravity could be very different 
at these scales from what we experience every day. It could be that the gravitational 
field is very slightly massive, and the fact was simply overlooked, because the effects 
are naively negligible at scales shorter than the characteristic Compton wavelength. 
But, it was shown [1J that a simple mass term, quadratic in the fluctuation of the 
metric, is not acceptable. The classical Hamiltonian - constructed in the ADM 
formalism - is not bounded from below. The model is not classically stable. 

Modifying gravity in the IR regime seems to be rather non-trivially constrained. 
In the past, several attempts to build classically stable models of IR modified gravity 
have been made. For instance, extra-dimensions may manifest themselves only 
at very large distances, like in the DGP model [2]. Or, without invoking extra- 
dimensions, the presence of non-local interactions could change the gravitational 
field in the IR regime, as in [3]. Or again, a similar effect can be induced by 
Lorentz- violating terms, like in [U El El El • 

In [8], we considered a class of models of the latter group. We were able to 
construct models that are four- dimensional and local, and have the property of being 
both stable in the IR and of departing from standard gravitational interaction only 
at large distances. But this was only possible at the price of explicitly breaking the 
Lorentz symmetry of the action, which had the effect of introducing new phenomena, 
like instantaneous interactions. These effects can be within the experimental bound 
by tuning the parameters, and more importantly they do not violate causality. 

In the present letter, we will extend the previous discussion by looking at ef- 
fective models over a constantly curved background. GR will be deformed with a 
small Lorentz-violating term. The model is classically stable and behaves like a 
de Sitter space-time at short scales. At larger and larger distances, the departure 
from de Sitter becomes more and more pronounced. 

A rather generic feature of the model is the tendency to stop the cosmological 
evolution and to lead towards a contracting phase. The effect of the deformation 
is opposite to that of the cosmological constant, for, as a cosmological constant 
increases the expansion speed, the deformation decreases it. After the expansion is 
stopped, i.e., when the Hubble parameter goes to zero, the Universe will go through 
a contracting phase, leading to a late-time (contracting) de Sitter phase. 

This picture will be modified by quantum corrections. In particular, out of 
the gravitational field, particles are produced when the Hubble parameter changes. 
The effect of quantum produced particles is to create at late times a space-time 
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singularity of the same kind as the one of a supercritical Universe. In the latter 
case, if the matter energy density is greater than a critical value, the scale factor 
grows up to a maximum size and then contracts towards a "Big Crunch" singularity: 
at some finite moment in the future the scalar curvature diverges. In the present 
case, instead, no criticality condition is present. The cosmo logical evolution is not 
stopped by the matter density, but by the deformation we introduced, therefore, no 
matter how few particles are produced by quantum effects, a singularity will always 
be reached. 

The model, which is classically stable, is unstable under quantum correction. 
This instability is milder, in fact its time scale can be made parametrically much 
greater than the present age of our Universe. 

The present letter is organized as follows: we will firstly discuss the general prop- 
erties of the gravitational field in models, where a Lorentz-violating interaction for 
the metric is added, explicitly showing the absence of additional degrees of freedom 
and the presence of instantaneous interaction. 

We will then study the cosmological solution of the (deformed) Einstein equation 
for two specific choices of the Lorentz-violating interaction. One will give rise to a 
bouncing Universe, while the other to a cyclic Universe of parametrically long period. 
In both cases, we will consider the quantum production that occurs at the bounce, 
and its effect on the cosmological evolution. 

The cyclic model will then be discussed in the context of inflation. We will find 
that the parameters of the model are constrained by the usual cosmological bounds. 
In particular, they have to be exponentially smaller than the Hubble scale during 
inflation for matching density perturbation with the values observed in the CMB. 

Finally, we will see the effects of the present class of deformations when in con- 
junction with the massive term studied in our previous work [8]. 

2 General Overview 

From a particle physicist's point of view, gravity is the interaction that emerges 
upon gauging the Lorentz symmetry of Special Relativity. As for any gauge theory, 
it should be possible to describe the low energy effective theory, which arises from 
the breaking of a part (or the whole) of the gauge symmetry as a result of some high 
energy dynamics. 

In gravity though, unexpected constraints arise for the low energy theory. It was 
noted by Boulware and Deser [1] that, have the classical non-linearities of the gravi- 
tational self-interaction taken into account, the Hamiltonian would generally not be 
bounded in models of massive gravity. Thus the Boulware-Deser (BD) instability 
appears. 

The origin of the instability - and therefore how to render stable a model of 
modified gravity as originally described in [8] can be easily understood in the 



3 



ADM formalism [9], that is the Hamiltonian formalism for gravity. 

Let us construct the GR Hamiltonian. By foliating the space-time with hyper- 
surfaces Ht f° r a time variable t, we can replace the four- dimensional metric with 
the following three-dimensional variables 

7^ = ^, N = (-(VT 1/2 and (2-1) 

N is known as the lapse function, and iVj as the shift function; 7^ is the induced 
metric on S t . In term of these variables, we can write the four-dimensional ones as 

= AV7, (2.2) 
Wft ee (3) 7£ + K,ijJC ij - K? , (2.3) 

where /Qj is the extrinsic curvature on S f , defined as 

*V = ^ [7« - ViJVj - V^i] • (2.4) 
The canonical momentum tt*- 7 ee bLjb^ij is related to /C^ by the relation 

7r« = ^ - £7"'] . (2.5) 

We now have all the ingredients to write the Hamiltonian for the Einstein-Hilbert 
Lagrangian 

C = y/gK -> W ee 7r%, - C\ iiS ^ = ^7 [A^° + , (2.6) 

where 

= - {3) n + 7- 1 (W j - ~tt 2 ) , 

i? ee -2Dj(-f- 1/2 n ij ) , (2.7) 
and Z)j is the covariant derivative defined with respect to 7^. 

Both iV and iVj appear linearly in the Hamiltonian, thus they are Lagrange multi- 
pliers. The variation with respect to them leads to the constraints - R° — and 
R l = - on the propagating degrees of freedom. The Hamiltonian is exactly zero 
on the surface of the constraints, hence the theory is trivially stable: the energy 
density of the system is bounded from below. 

We will show the emergence of the BD instability for deformed Einstein-Hilbert 
actions. For sake of definiteness, we shall consider the Pauli-Fierz (PF) model [TO] . 

The PF term is the most generic deformation that is quadratic in the fluctuation 
of the metric over a particular background, and it is Lorentz-symmetric. It describes 
a mass for the gravitational field 

- \mlAK v ' (^) 2 ] = -\ml F [h% ~h 2 - 2N 2 + 2h(l - N 2 - N 2 )] . (2.8) 
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In the equality, the field h^ v is expressed in terms of the three-dimensional variables 
( 12. ip . The tensor field is defined as h^ u = g^ v — g^ u over a particular background 
metric g^ v . The indices are contracted using g^ u , i.e., = g^h^, h = g^hij, and 
so on. 

It is evident that the lapse and the shift functions cease to be Lagrange multi- 
pliers, and the variations with respect to them lead to algebraic equations for them, 
rather than constraints on the propagating degrees of freedom, as in GR. This is 
hardly unexpected. A massive field is known, from the Lorentz group representa- 
tion, to propagate a number of degrees of freedom different from that of a massless 
field. 

But, if we now study closely the equations 

2rrip F h ' ( ) 

N l = -l-^i-^)- 1 ^, (2.10) 

we notice that the Hamiltonian, after substituting the above values for N and iVj, 
h 



H 1 



4m| F 



+ l m ( n ,(h; i -lr + 2h). (2.11! 



is unbounded, as it is readily seen by considering the limit h — *■ , while keeping 
^fR° and R i = fixed. 

It appears that models with iV 2 -terms - such as PF - are generally unstable. 
An easy way-out is obviously to consider more general classes of deformations, in 
particular the ones linear in N. It should be noted though, that the iV 2 -term in the 
PF Hamiltonian (12. 8p comes from the time component of the tensor h^. From the 
definition of the lapse function (12.11) . it follows that h 00 = g 00 — g 00 ~ iV 2 . Thus, 
removing such a term would lead to an explicit breaking of the Lorentz symmetry 
in the action. For a detailed discussion of a PF-like model with such a property we 
remind to our previous letter [8]. 

The class of models we would like to discuss in the present work has the following 
Hamiltonian 

U = [NR° + NiR 1 + 2AiV - 2m 2 iV/^)] , (2.12) 

where /(-v/7) is some function of the determinant of the spatial metric 7^ as in 
( 12. ip . and we have assumed the presence of a cosmological constant A. The defor- 
mation we added is modelled to be linear in N, thus the lapse function is still a 
Lagrange multiplier. As in GR, the Hamiltonian is exactly zero on the solution of 
the constraints. The model is therefore (classically) stable. 

This model, as well as the ones presented in [8], should be thought of as an 
effective low energy theory. Like in the Higgs mechanism, the gauge symmetry is 
broken at low energy, hence we are assuming the presence of some UV-physics that 
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spontaneously breaks the Lorentz symmetry of the action, like in the recent models 

euni-ini- * 

The introduction of a deformation could, in principle, lead to the propagation 
of more degrees of freedom, some of which may develop into instabilities for the 
theory. We will show that only a tranverse-traceless tensor mode is propagating. 
The explicit breaking of the Lorentz symmetry will instead show up as instanteneous 
interactions, as first noticed in To explicitly study the degrees of freedom, we 
turn now to the Lagrangian formalism. 

The Lagrangian can be found by performing a Legendre transformation on (I2.12p 

C = [U - 2 A + 2m 2 /( v / 7)] , (2-13) 

and it should be noted that we are forced to keep a somewhat mixed formalism. The 
deformation is written in terms of the determinant of 7y, hence retaining in part 
the notion of the three-dimensional variables (12.11) used in the ADM formalism. 
The equations of motion are 



^-m 2 ^y ) ^ = 0, (2.14) 



where G^ u is the Einstein tensor defined as G^ u = 1Z^ V — l/2TZg^ u , and the last 
term is zero for //, v ^ 0. 

The deformed action ( 12.131) is not invariant under Lorentz transformations any- 
more. The determinant of 7^ will transform under — > x a as it can be explicitly 
checked. The introduced deformation breaks the Lorentz invariance of the action 
down to the rotational group. The measure on the hypersurface S t is invariant 
under diffeomorphisms acting on its own world-volume. 

The breaking of the Lorentz symmetry stems out from the presence of a preferred 
frame in the model. In the construction of the GR Hamiltonian, a frame is chosen 
when picking up a particular foliation Ej of the space-time. The choice is formal 
in GR, where the action is invariant under the entire group of the diffeomorphisms, 
but it is not in the theory at hand. 

The deformation we consider is a function of the spatial part 7^- of the metric 
tensor only. To define 7^, we need to choose a direction and 7^ be the induced 
three-dimensional metric on the hypersurface orthogonal to u^. The preferred frame 
is defined therefore by u^. Any transformation orthogonal to is a symmetry of 
the model, as briefly stated before. 

Because of this remaining invariance of the action, some components of the 
perturbation over a background metric g^ u are not independent to the others. 
For sake of simplicity, we shall consider g^ v to be the Minkowski metric rj^. It can 
be shown that the results we are presenting do not depend on this choice. 

The redundancy in the components of h^ u can be fixed by a gauge choice, for 
instance d % hij = 0. It fixes the longitudinal modes of hij, but not its trace. The 
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action is indeed not invariant under the transformation h\ — > h l i + 2<9*£j, hence the 
trace h\ can not be fixed by a gauge choiceEl 

The gauge fixing is not the only condition we have to impose on h^. The 
requirement that the equations of motion are covariantly conserved 

(2f + V7/'UV7 + V a d a ^ - 9 a "(r%5l + r^°)^ = , (2.15) 

which is known as Proca condition for massive vector fields, constraints even further 
the independent components of h^ v . 

For perturbations over a Minkowski background, ( 12.151) becomes 

&d i h Qi = 0, (2.16) 
(/o + ^m^ + ^o^°o = 0, (2.17) 

where = /'(Vt) and f» = f'(Vl)- 

The study can be done easily using the following decomposition, common in the 
study of cosmological perturbations [19] . 



v % B vvx + difyE j ' ' v v ^ ; ' Miy V ^S T 

where /i^ T is a transverse and traceless tensor, ipi, Fj transverse vectors, and the 
rest scalars. 

The scalar E and the vector Fj are fixed by our choice of gauge d l hij = 0. 
The analogous (12.161) . (12. 17ft of the Proca conditions fix instead the longitudinal 
component of hoi, i.e., the scalar B, and a combination of h° and h\, that is the 
solution of (I2T?D . 

The independent components of h^ v are therefore: a transverse-traceless tensor 
hJ 3 T , a transverse vector ipi, and a scalar, combination of <ft and 

Not all of these five independent components are propagating degrees of freedom. 
It was noticed in [S] that a common feature of Lorentz-violating theories is the 
presence of an instantaneous interaction. An explicit study of the equations of 
motion for each and all components of the perturbation would show that both 
the vector and the scalar appear without time derivatives in their equations of 
motion. Hence, they can not be identified as propagating degrees of freedom, rather 
as an instantaneous background. 

The propagating degrees of freedom of the model are equivalent to the ones of 
standard GR, i.e., a transverse-traceless tensor field. The difference is the presence 
of instantaneous interactions. It should not surprise, in fact breaking the gauge 



*Alvarez et al. in [18] studied a similar model, where the action is symmetric under transverse 
diffcomorphism (TDiff), that is h^ u — * h^ + d^^) with <9 M £ M = 0. They found that TDiff invariant 
theories contain an additional scalar field; in our case the action is invariant under spatial TDiff 
and no additional propagating degrees of freedom are present. 
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invariance of the action does not allow to remove from scattering cross sections 
components of the metric, like the Newtonian potential h° , which is instantaneous 
in nature. Thus, they would appear not only in the exchange of virtual particles, 
but they would also manifest as physical phenomena. 



3 Stopping the Cosmological Expansion 

Having established the properties of the model under discussion, we shall now con- 
sider more specific examples. We choose /(-v/7) = 7 a//2 . 

We are interested in particular to IR-modifications of gravity. We would like the 
deformation to be dominant only at late times during the cosmological evolution, 
i.e., when a(t) ^> 1. This can be achieved by assuming a > 0. 

The equations of motion for this particular choice of fiy/^y) are 

Gfw + [A - (a + 1) m 2 7 Q / 2 ] g, u - a m 2 5° Sjj = . (3.1) 

On the FRW ansatz with zero spatial curvature, ds 2 = — di 2 + a(t) 2 dx 2 , they 
become 

d 2 A m 2 ,„ , 

the space-space component of (13 .ip is proportional to the time derivative of (13 .2p as 
it can be checked explicitly. 

The additional term a 3a in the Friedman equation can be mimicked by a field 
with equation of state w = —1 — a. For a > 0, this fictitious field has w < —1. 

When m 2 > 0, the energy density of the "field" has a wrong sign, while if m? < 
its energy density is well-behaved. The former case will be discussed in detail in 
the following of the present section. The latter instead, being of a matter field with 
w < — 1, will drive the cosmological expansion at an ever increasing acceleration rate, 
as opposed to de Sitter space-time of constant acceleration rate. Such a matter field 
(or deformation) will drive the Universe towards a final state sometimes called "Big 
Rip" : the scale factor will diverge in a finite time. 

Our intuition on the cosmological evolution for the m 2 > case is scarcer, and 
this case should not be treated as for a matter field. Let us therefore investigate the 
solution step by stepQ. 

It might be proven useful to discuss the solution of the (modified) Friedmann 
equation (13. 2p . as of a classical point-particle moving in a potential V(a) = (— A/3 + 
m 2 /3a 3a )a 2 with zero total energy. The potential is pictured in Fig.fH and it should 
be noticed the presence of a turning point at a* = (A/m 2 ) 1 / 3 ". 

A classical point-particle moves down the potential hill from an initial position at 
a = until stops at a*, because of the "attractive force" generated by the m 2 -term. 



T The study of the cosmology dynamics in the presence of a matter field with negative energy 
density can be found in [20] . 
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V(a) 




Figure 1: Potential for the analogous 1-dimensional dynamics, the dashed line is the 
potential for m — 0. 

After that position is reached, it rolls down in opposite direction towards a = 0, 
which is reached in an infinite time. 

This classical analog is easily translated into the cosmological evolution of the 
Universe. For small scale factor a ~ 0, the dynamics is dominated by the cosmo- 
logical constant and the Universe is in an approximate de Sitter phase. During 
this time of exponential expansion, the "attractive force" of the m 2 -term will grow 
in intensity, until it will become dominant driving the cosmological evolution to a 
bounce at a(t*) = a* = (A/m 2 ) 1 / 3 ". Then, it will contract approaching at late times 
a (contracting) de Sitter phase. 

We can solve analytically the equation of motion (13.21) 

2 

where to is an integration constant to be fixed by imposing the initial condition 

a\t=o = 1: to = -j^t- cosh -1 */ ■Ar. 

This exact solution is plotted for some positive value of a alongside with the Ricci 
curvature TZ(t) in Fig. [2j The features previously described are easily recognizable: 
at early and late time, the expansion is dominated by the cosmological constant and 
the Universe exponentially expands and contracts respectively. 

The dynamics presented here is classical, and would be modified by quantum 
corrections. In particular it is known that, whenever the scalar curvature changes 
in time, particles are created via a phenomenon similar to the Hawking radiation 
for a Black Hole. 

To understand it, let us consider a basis of particle creation/annihilation oper- 
ators {v4|,,y4fc}. The Fock space is defined by determining the vacuum state \fl), 
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a(t) 



R(t) 




t 



Figure 2: Plot of a(t) and K(t). 



which is destroyed by all the annihilation operators, Afc|f2) = 0, and then populating 
it by acting with the creation operators onto 

The operators are constructed by canonically quantizing the fields of the par- 
ticular theory at hand, but in doing so we should assume a particular background 
metric. Hence, if the metric is time-dependent, so are the operators. The vacuum 
\Q) is annihilated by all Ak at a given time, but not, in general, at every time. 
Thus as time goes by, the vacuum state will be in a superposition of particles, for 
A k (t)\Q) / at a generic late time. 

Since the background we found is time-dependent, we would like to study the 
effects quantum particle production has on the cosmological evolution. 

For sake of clarity we specialize to a = 1. Following |21j, the energy density p q 
created up to a time t is given by the following expression 



where V'(t) = (1 — 6£)(d 2 + ad) for the FRW ansatz and £ = 1/6 for conformally 
coupled fields; rj is the conformal time defined as dr] = dt/a(t) 



and I(z;a,b) is the regularized beta function I(z; a, b) = B(z;a,b)/B(a,b), with 
B(z; a, b) and B(a, b) being the incomplete and complete beta function respectively. 

For sake of simplicity we approximate the quantum corrected dynamics as if 
particles were generated all at once when TZ(t) = 0. This rough approximation will 
not change the qualitative picture we will describe. 

By evaluating the integral of (13.41) up to £, we find the energy density to be 




(3.4) 




(3.5) 



1 



A 2 ~ A 2 , 



(3.6) 



288tt 2 



where X is the numerical result of the (adimensional) integral of ( 13. 4ft . 
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a(t) 




t 



Figure 3: Evolution of the scale factor considering particle creation at t. 



As it was noticed in [21] . p q does not depend on the time when the de Sitter 
expansion stops, but solely on the change in the scalar curvature. In our model, this 
means p q is proportional to A, but not to m 2 . 

After t, the cosmological evolution will change to include the presence of matter. 
The analytical solution of the Friedmann equation is 

- [™ V v m^) [f^^- i-yrr^J ■ (3 ' 7) 

where sn(w, //) is one of the Jacobi elliptic functions. 

The plot of the scale factor a(t) obtained by matching (13. 3p with ([3.7]) at t — i 
is shown in Fig. [3] 

At t > t , the m 2 -term is driving the cosmological evolution, in fact to have it 
to stop the exponential expansion driven by A - that is to approach TZ(t) = - 
its "strength" has to be of the same order of the cosmological constant. Because 
m 2 is dominant, the cosmological evolution does not change significantly from the 
one we described without matter. The scale factor reaches a maximum, and then 
starts to contract. The dynamics starts to depart from (13. 3p at this point. Instead 
of approaching a contracting de Sitter phase with constant curvature A, at a distant 
but finite time in the future the Ricci curvature diverges: a "Big Crunch" occurs. 
At late times, when - during the contracting phase - a(t) <C 1, the matter density 
is driving the evolution, and therefore the singularity is unavoidable. 

This picture is similar to a supercritical Universe, with the difference that no 
criticality condition for p q is present. The expansion is not stopped by the (super- 
critical) matter density, but by the IR modification we introduced. Thus no matter 
how small p q is, the dynamics ends with a singularity at a finite time. 
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The model, despite being classically stable, is unstable under quantum correc- 
tions. 

Moreover, even considering that the "Big Crunch" will occur in a (paramet- 
rically) very distant future, the exit of the de Sitter expansion, which could be 
associated with a period of inflation, coincides with an era dominated by the m 2 - 
term. Thus, the model would not follow the known cosmological evolution of our 
Universe. 



4 Cyclic Universe 

As we have discussed in the previous section, a "Big Crunch" is unavoidable once 
quantum corrections are taken into account. One way to avoid the singularity could 
be to modify once more the gravitational Lagrangian. 

By studying the effect the m 2 -term has on the cosmological evolution, we under- 
stood that the modification we introduced slows down the expansion and ultimately 
stops it, after which the Universe goes through a contracting phase. If we could 
have a term that acts as the m 2 -term, namely that slows down the cosmological 
evolution, but is dominant only when the size of the Universe is small, we might 
stop the collapse before the scale factor reaches zero, and the curvature diverges. 
This can be attained by having chosen the exponent a to be negative. 

The most general (and minimal) Lagrangian is 

C = [R, - 2A + 2mV /2 + 2fcV /2 ] , (4.1) 

where a is positive and (3 negative definite. On the FRW ansatz, the equation of 
motion is 

a 2 m 2 o„ k 2 o„ A . 
- + —a 3a + -a 3 ?-- = 0. 4.2 

GT 3 3 3 

We can study the dynamics of a(t) by analogy with a classical particle moving 
in the potential V = (m 2 a 3a — A + k 2 a 3l3 )a 2 /3. The potential is plotted in Fig. BJ 
for (3 < —2/3, (3 = —2/3 and (3 > —2/3. The common feature of all the plotted 
potentials is the presence of two turning points a(±), between which the point-particle 
would oscillate back and forth. 

The model is of a cyclic Universe that "eternally" oscillates between a minimum 
and a maximum size. Between the two extrema, the Universe is in a de Sitter phase, 
either contracting or expanding. 

We can study analytically the dynamics for the choice of a = 2/3 and (3 = —2/3; 
the solution, shown in Fig. is 



a(t) = or-) 




' 2 2 
a (-) ~ 



1/2 



(4.3) 
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Figure 4: Potential for the analogous 1-dimensional dynamics; the continuous line 
is for p = —2/3, the dashed for (3 > —2/3 and the dot-dashed for (3 < —2/3. 



where ac±\ = (A±VA 2 — Ak 2 m 2 ) /2m 2 , and sn(w, li) is the Jacobi sn elliptic function 
with periodicity!*! 



4v/3 

T 



(A 2 - Ak 2 m 2 y/ 4 




+ K 



a l+) 



a 2 - a 2 

"(-) "(+) 



(4.4) 



where is the complete elliptic integral of first kind. 

The presence of the fc 2 -term effectively "screens" regions of small size. It slows 
and stops the contraction up to a non-zero scale factor O(-), in the same way that 
the m 2 -term stops the expansion at large scale. 

Naively we would expect this model to be stable under quantum correction, for 
/c 2 -term could stop the collapse in presence of matter too. It is like the new term 
creates a potential barrier at small scales, preventing the scale factor to reach zero 
size. The presence of matter would lower the barrier, but by tuning k 2 we can avoid 
its disappearance. 

The loop-hole in this argument is the periodicity of the model. At each cycle 
new matter is generated by quantum corrections as seen in the previous section. 
Eventually p q will "overcome" the barrier leading towards a "Big Crunch", unless 
the /c 2 -deformation is always dominant at small scale factors, no matter how big p q 
is. This is obtained by taking (3 < — 1 — w, where w is the equation of state for the 
particles generated through quantum effects. For radiation w is equal to 1/3, thus 
if /3 < —4/3 a minimal size 0(_) is always present. 



tmore exactly, the Jacobi sn(u,/i) is doubly periodic in the complex plane, that is sn(u + 
4(if(/i) + i n K(l — fj,)), ji) = sn(u, fi) where n£ Q. In the present case, n is fixed by requiring the 
period to be real. 
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a[t] 




a 



Figure 5: Plot of the scale factor a(t) in respect to time for the cyclic solution. 

In the approximation the m 2 -deformation and cosmological constant are negligi- 
ble, which is always the case for a(t) <C 1, the minimal size is 



Obviously, because more and more matter is generated at each cycle, a time will 
come when 0(_) will be of Planck size, and therefore our semi-classical description 
will break down. 

5 Embedding Inflation 

In the following, we would like to depart slightly from the previous discussions, that 
have considered the effects of gravity modifications in a de Sitter space-time. 

We will consider a generic inflationary model, and describe how the terms we 
introduced would effect inflation, and what kind of bounds we could have. We will, 
therefore, consider modified gravity coupled to a scalar field, the inflaton. The plot 
of a typical inflaton potential V((p) is sketched in Fig. [61 

When the inflaton is atop the plateau, its energy density is dominated by the 
potential, namely ip 2 <C V(ip), and approximately constant. This plays the role of 
the cosmological constant A = V(ip). The space-time is exponentially expanding as 
long as the inflaton sits on the plateau. 

This configuration is not stable, in fact it is energetically favorable for the scalar 
field to condense at the minimum of the potential, (<p) = <pq. When the inflaton 
condenses on its true vacuum, inflation ends, because the energy density is no longer 
a non-zero constant. 




(4.5) 
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V[ip] 




Figure 6: Typical inflaton potential V(<p). 

Once inflation is over, the inflaton decays more or less efficiently, depending on 
its coupling, into the Standard Model fields, reheating the Universe. From this point 
on, the cosmological evolution follows the Standard Model of Cosmology 

We can embed this general picture into our model and describe what constraints 
arise from cosmology and general requirements. 

The Universe described by (14. 3 p - solution of the Friedmann equation - goes 
through periodic exponential expansions and contractions. The requirement to have 
the usual description from the inflation framework within this model demands the 
scalar field <p to condense into its true vacuum during an exponential expansion. 
Conversely, if it condensed during a contraction, the density perturbation originated 
during the inflationary period would be washed out at the unavoidable bounce. 

Thus the inflaton can condense any time between a(_) and ci(+), and this in turn 
provides a bound on the number of e-foldings N e as a function of the parameters 
m 2 and k 2 

G (+) i A + VA 2 - 4fc 2 m 2 

N e > log = log (5.1) 

O(-) A - VA 2 - 4Pm 2 

This relation can also be read as a bound on the parameter m 2 . Assuming m 2 
to be much smaller than A and k 2 , we find the following bound 

m 2 <^exp[-2iV e ]. (5.2) 

The modification we introduced must be exponentially smaller than the de Sit- 
ter curvature in order to satisfy the cosmological bound on N e as it arises from 
matching the inflaton density perturbations with the anisotropics observed in the 
CMB (iV e >60). 
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Few comments should follow. Firstly, (15.21) is strictly an upper bound, for, if 
inflation end when the m 2 -term is dominant, that is when a(t) ~ «(+), the cosmo- 
logical evolution could not start in a radiation-dominated epoch, as in the Standard 
Model of Cosmology. 

Secondly, whether in our model inflation could be eternal - see, for instance 
[22| ,23j. As we have throughly discussed in the previous sections, the additional 
potential-like term acts to stop the cosmological evolution bringing to an end any 
inflationary periods. It is not obvious how inflation could be eternal in a model that 
does not allow the space time to expand indefinitely. 

In the standard picture of eternal inflation, false vacuum bubbles nucleate due 
to quantum fluctuations of the inflaton field. If the size of these bubbles is larger 
than their Hubble radius, they are causally disconnected from the "ambient" space. 
They will independently evolve and eventually pinch off. 

As long as the deformation is negligible, the previous picture applies to our model 
as well. From the point of view of an observer living in one of these bubbles, they 
have no way to know anything about what is happening outside the horizon. It 
should be natural to expect that, once the "ambient" space bounces back, it keeps 
collapsing towards a de Sitter contracting phase, while the bubbles, unaware of 
anything happening outside their Hubble radius, expand until they pinch off. Hence 
new "baby" Universes are generated, and they will follow the same evolution of their 
"parent" Universe: the Universe will keep self-replicating and at any time at least 
one patch will be in an inflationary regime. 

Even though every bubble can expand only till reaching a finite size before bounc- 
ing back towards a contracting phase, the total number of nucleated bubbles will be 
infinite, hence eternal inflation is a possible scenario in our model. 

5.1 A Curiosity 

Something curious happens for a particular, non-small value of km/ A. 

If A = 2km, the two turning points a(±) are equal as it follows from their defi- 
nitions. At this particular value, the solution of the Friedmann equation is of a flat 
space-time. 

For the tuned value of m 2 , an otherwise de Sitter space-time turns out to be 
effectively flat. This solution can be nicely understood from the analogous classical 
point-particle description. As we have already stressed, the cosmological constant 
acts as a repulsive force, whilst the m 2 -term as an attractive one. When A = 2km, 
those two forces balance exactly leading to the allowed range of a(t) to shrink down 
to a point. From the cosmological view point, the scale factor is time-independent 
and therefore the space-time is effectively flat. 

Some fine-tuning is required for this solution, therefore it would be interesting 
to understand how (if) the solution is reached in a generic inflationary model. 

Let us consider a linear potential V(<p) = A — rfip with rf <C A to satisfy 
the slow-roll condition. In a usual inflationary model, inflation would be eternal, 
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for no minimum of the potential of the scalar field is present. But an exponential 
expansion always stops when the m 2 -term becomes dominant, so it is not clear what 
dynamics will follow in our model. 

The scalar field follows its own equations of motion coupled to gravity 



(p + 3-cp + V'(cp) = 
a 



a i 2 2 , ; 2 -2 

— + ma + k fa — 
a 2 



\v 2 + V{ V ) 



0, (5.3) 



where H = a/a, the Hubble constant, is a friction force for ip. In the slow-roll 
approximation, the inflaton reaches critical but small velocity \<p\ ~ \V'((p)/3H\ = 
r) 2 /3H <£V(<p). 

The effective cosmological constant is A ~ V(<p) = Ao — i] 2 Lp. Since <p is changing 
in time, so is A, but we can assume the slow-roll approximation to hold at any time, 

The Hubble constant H 2 = (a/ a) 2 is monotonically decreasing from its initial 
value Ao- The lower limit is at A ~ 2km, where (p diverges. Thus, at this point, the 
inflaton is in a fast-roll regime with the kinetic energy dominating over the potential: 

We can drop the potential from the equations of motion (15.31) 

<p + 3-ip = 0, (5.4) 
a 

a 2 1 

— + mV + k 2 a- 2 - -v? 2 = 0. (5.5) 
a z 2 

The dominant term is the kinetic energy - ip> 2 = <p>la~ 6 from H5.4j) - that drives 
the expansion until, at very late times, the m 2 -term becomes dominant and leads 
to a contracting phase. 

Because (p\ ^> m 2 , k 2 , the dynamics is that of a supercritical system, and the 
Universe would eventually collapse in a "Big Crunch" singularity. 



6 Massive Modification 

In the present section, we will describe the effects a massive modification - of the kind 
first described in [S] - has when it is considered in conjunction with the deformation 
described in sectional 

Before starting on describing the cosmological evolution that arises when both 
deformations are taken into account, we should emphasize that the present modifi- 
cation is very different from the one so far described. We showed that the number of 
propagating degrees of freedom is unaltered when sJ—gji^Ff) is introduced. Instead, 
in the case we shall present the number of degrees of freedom will be different. 
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The Lagrangian of the model we would like to discuss is 



C = y/=g[K - 2A + 2/4 
with the following equations of motion 



N 



+ 2/47 



(6.r 



G 



A 



2 1 



N 
-I) 



a/2 



9n» + 
= 0. 



(6.2) 



Notation is consistent with the one used before. N is the lapse function and 7^ the 
induced spatial metric defined in (12.11) . As we discussed in our previous letter [8] , the 
modification y/— g(l — N) 2 /N is the only one quadratic in the lapse function, over 
which fluctuations do not present instabilities like tadpoles and the Hamiltonian of 
(16. ip is bounded from below. 

It is evident from the Hamiltonian 



H = y/j [NR° + NjR j + 2AN - 2/i 2 (l - A^) 2 - 2^ 2 2 N 1 a/2 ] 



(6.3) 



that the shift function Nj remains a Lagrange multiplier, while the lapse A" ceases 
to be it. The algebraic equation for N, that is 2[i\N = R°/2 + A + ^\ - ^l^ a/2 
can be seen as a constraint for the degrees of freedom, thus on this background the 
gravitational field propagates three degrees of freedom, instead of the two as in the 
previous cases. 

For this, the newly added modification is different from the ones we discussed 
previously: despite our analysis will follow closely the one of the previous sections, 
we should bear in mind that the models describe two very different fields. On the 
one hand, in section [3] the field is a massless tensor field propagating on a deformed 
(with respect to standard GR) background; on the other, in the present section we 
will discuss a model for a massive^] tensor field. The presence of the deformations 
will modify the cosmological evolution away from GR, as we will readily see in a 
moment. 

On the ansatz ds 2 = — A^(t) 2 dt 2 + a(t)dx 2 and after a bit of algebraic manipu- 
lation of the equations of motion (I6.2p . we find the following equation for the scale 
factor a 

.0 A 



1 + 2^(1- Vl + a" 3 ) 



A 



— a 3a 
A a 



0. 



(6.4) 



where the derivative a is in respect to the proper time dr = N dt and Af 
relation that can be derived directly from the equations of motion. 



§in this context, we call "massive" a field that propagates a number of degrees of freedom 
different than two; it should be noted also that the action is not symmetric under the Lorentz 
symmetry and therefore the intuition of a massive tensor field propagating five degrees of freedom 
is not necessarily respected. 
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V[a] 




Figure 7: Potential for the analogous 1- dimensional dynamics: continuous line is for 
A*i = A*2 = 0, long dashes for /i 2 = 0, short for /xi = and mixed dashed line for the 
general case. 



It is useful to discuss the solution of H6.4[) as for an analogous one-dimensional 
point particle moving in the potential V(a) = — fo 2 1 + 2^(1 — VX + a~ 3 ) — ^-a 3a 

The potential is depicted in Fig. [7] for various values of /^i and \i2- In particular, 
the continuous line is for [L\ = ^ — and the dynamics is that of a de Sitter 
space-time with cosmological constant A, as expected. The other two limiting cases 
are when one of the two deformations is absent; we find that for \L\ = - short 
dash curve - the cosmological expansion proceeds as for a de Sitter universe until a 
maximum size 0(+), after which it goes through an exponentially contracting phase. 
That is, the solution we discussed previously in so much detail is recovered. 

For /i 2 = - long dash curve - a minimum size a(_) for the scale factor emerges. 
The dynamics is very much similar to the one we described for the cyclic universe 
model of section H] when only the fc 2 -deformation of (14. Xp is present: the cosmological 
evolution is that of a contracting de Sitter universe until the Universe reaches size 
at-), after which it bounces back toward an expanding de Sitter phase. 

When, instead, both deformations are present, the analogous potential V(a) 
displays two turning points am. The dynamics can be read easily from it, and 
it is like the cyclic Universe model of section HJ the cosmological evolution goes 
through subsequent periods of expansion and contraction, between the minimum 
and maximum sizes ar_) and 0(+). Those values are the zeros of V(a), and given 
that, when the size of the Universe is of order of a(_) the /i|-term is sub-dominant 
in respect to both the cosmological constant A and the /i 2 -term (and vice-versa for 
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a(+)), they are approximately: 

( *j4 V 

/ A \ ^ 

a (+) ~ (^j . (6.5) 

The turning point a(+) is the same as the one we found in the model of section [31 
On the other hand, ~ (fix/ A) 2 / 3 , when /i 2 -C A, is what would be the turning 
point for a deformation of the kind j a ^ 2 with a = —1/2. This is not surprising. 
In the regime when /i 2 is sub-dominant, the cosmic evolution is driven by the fif- 
deformation. From the equations of motion f!6.4p . we find cosmology being driven 
by a term - \/T+ a -3 ~ a -3 / 2 - that acts like a 3a for a = —1/2, i.e., equivalent to 
the deformation of section [3] for a particular choice of a. 

Again we should stress that, even though the backgrounds are alike, the fields 
propagating over them act very differently, for one carries two degrees of freedom, 
while the other three. Therefore, they describe two very different gravitational 
models. 



7 Conclusions 

Models, that describe gravity beyond Einstein GR, have been focussed mainly on 
higher derivative deformations, like f(TZ) gravity for instance. In the present letter, 
we discussed a class of models, that modify gravity via potential-like terms. 

The /(a/7) deformation we introduced does not add any derivatives of the metric 
to the action. Because of this, we argued that the propagating degrees of freedom 
are of a transverse-traceless tensor field, as in GR. We arrived to this conclusion by 
studying the perturbations of the metric in the Lagrangian formalism. Therein, it 
was also evidenced the presence of instantaneous interactions. This a-casual effect is 
characteristic of the models at hand, where Lorentz symmetry is broken explicitly. 
Their presence was firstly noted in [8], and we remind to it for a more detailed 
discussion. 

We then studied the exact cosmological solution for some particular choices of 
/(•v/t). The main feature is that the introduced term acts to generically stop the 
cosmological evolution. Depending on the details of the model, a bounce is generated 
during either a period of contraction or one of expansion. Independently from these 
details though, particles are produced via quantum effects at the bounce. Quantum 
corrections destabilize an otherwise stable classical solution. The result is to create a 
future "Big Crunch" singularity: after a finite time the scalar curvature will diverge 
and the Universe will shrink to zero size. 

In more specific sections, we discussed the effects and bounds on the parameters 
that arise when the model is embedded into a generic model of inflation. In par- 
ticular, we noticed that the dimension-full parameter m 2 has to be exponentially 
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smaller than the Hubble radius during inflation, so to satisfy the bounds imposed 
by the CMB. 

Left to future investigations is the study of Schwarzschild solutions for the present 
model, like in [2U (25J [26] . It would be interesting to know what kind of effects the 
/(y / 7)-term has on a Schwarzschild-like solution, and to see how (if) it screens the 
gravitational field of a massive point-particle. 

Also left out is the understanding of the UV completion of the model. We 
have always stressed that the introduced deformation /(v/7) should be thought 
as an effective term rising from some UV phenomena. The question is to exactly 
determine what kind of phenomena. Studies in this direction have flourished in the 
past few years, [I] and [H]-[I7], and it would be interesting to have a fully consistent 
mechanism that breaks spontaneously the Lorentz symmetry and generates f{^)- 
term at low energies. 
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